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We discuss the possibility of extending the RG flows to complex coupling spaces. We argue 
that the Fisher's zeros are located at the boundary of the complex basin of attraction of IR fixed 
points. We support this picture with numerical calculations at finite volume for 2D 0(N) models 
in the large-rV limit and the hierarchical Ising model using the two-lattice matching method. We 
present numerical evidence supporting the idea that, as the volume increases, the Fisher's zeros of 
4-dimensional pure gauge SU(2) lattice gauge theory with a Wilson action, stabilize at a distance 
larger than 0.1 from the real axis in the complex j3 = A/g 2 plane. We show that when a positive 
adjoint term is added, the zeros get closer to the real axis. We compare the situation with the U (1 ) 
case. We discuss the implications of this new framework for proofs of confinement and searches 
for nontrivial IR fixed points in models beyond the standard model. 
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1. Introduction 

There has been a renewed interest in the Lattice community |2|, ]3|, Q, ||, |6| 0] for the possi- 
bility [Jsj] of having nontrivial infrared fixed points in asymptotically free gauge theories (see recent 
reviews by deGrand [[J and Ogilvie [[H]]). One particularly interesting situation from a phenome- 
logical point of view is when the j8 function approaches zero from below and encounters only 
small changes over a significant range of scale. We then say that the "running" coupling constant 
"walks". From this perspective, the ability to control the height of the j8 function appears to be 
quite important. A simple model where this can be done easily is the case of a quadratic j8 func- 
tion which appeals naturally in the quantum mechanical l/r 2 potential and other problems where 



conformality can be lost and restored [ ]1 1| , |12| ]. For a sufficiently large value of the constant term 
of a quadratic /3 function, one infrared (IR) and one ultraviolet (UV) fixed points are present. By 
lowering this constant term, we can make the two fixed points coalesce and then disappear in the 
complex plane. The situation is sketched in Fig. [j]. 



This motivated us [ |1 3[ ] to study extensions of renormalization group (RG) flows in the complex 
coupling plane. A general feature that we observed is that the Fisher's zeros - the zeros of the 
partition function in the complex coupling plane - apparently act as "gates" for the RG flows ending 
at the strongly coupled fixed point. This can be seen as a complex extension of the general picture 
proposed by Tomboulis [|14|] to prove confinement: the gate stays open as the volume increases and 
the flows starting in a complex neighborhood the UV fixed point can reach the IR fixed point where 
confinement and the existence of a mass gap are clearly present. In general, losing conformality 
corresponds to the generation of a mass gap and the presence of confinement and complex fixed 
points not on the real axis. We argue that such fixed points are related to the absence of Fisher's 
zeros on the real axis. 

In the following, we illustrate this scenario with model calculations for the 2D 0(N) non- 
linear sigma models in the large-A^ limit and the Ising hierarchical model. In all cases, the Fisher's 
zeros (of the partition function) seem to govern the global behavior of the flows near the real 
axis. In the infinite volume limit, these zeros delimit the boundary of the basin of attraction of 
the strongly coupled fixed point. For confining models, a "gate" remains open. We considered 
modifications or deformations that may affect that behavior (finite volume, change of dimension, 
additional pieces in the action). We will then present recent results regarding the Fisher zeros for 




Figure 1: Sketch of what happens when we reduce the constant term in a quadratic j3 function: the IR and 
UV fixed points merge and disappear in the complex plane 
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U(l ) and SU (2) 4D LGT (zeros at different volume; no RG flows yet). Note that at finite volume 
all the models considered here have a partition function analytical in the entire complex j8 °< 1/g 2 
plane. We should also mention previous studies of RG transformations in the complex temperature 



or coupling plane for ID spin models and 2D gauge models [|15fl. 



2. Complex RG flows in spin models 

We now discuss numerical calculations of complex RG flows in spin models. We first consider 
the 2D 0(N) non-linear sigma model in the large-Af limit. The partition function reads: 



(2.1) 



We use the notation /$ = 1/ (g^N), not to be confused with the j8 function, for the inverse 't Hooft 
coupling and M = m gap /Auv for the mass gap in cutoff units. For large-Af, we have the gap 
equation: 

^ - 1 1 d k 1 (2.2) 



^ M ' J- Jl J- 7 t(2K) 2 2(2-cos(k l )-cos(k 2 ))+M 2 ' 
At small coupling we have the asymptotically free relation: fi(M 2 ) ~ 1/ ln(l/M 2 ) 



Flows of 0(N) modal in infinite V 




Figure 2: RG flows obtained by rescaling complex values of M in Eq. (2.2) 



We first consider a simple complex extension of the running of the coupling when the UV 
cutoff is decreased: we take the mass gap on a small circle in the complex plane: m gap = ee' e and 
then lower the cutoff. A^y — > Aj/y / b. The flows for b = 2 are represented by red arrows in Fig. ^| . 
The blending blue crosses (called "blue lines" hereafter) are the j3 images obtained from Eq. ( |2~2| ) 
of two lines of points located very close above and below the [—8,0] cut of /3 (M 2 ) in the M 2 plane. 
The RG flows are constrained to stay inside the blue lines and we will explain why the Fisher's 
zeros are expected to stay outside the blue lines. The flatness of the flows at large j8 follows from 
the asymptotic freedom running: the complex phase in M 2 is not affected by the cutoff and its 
logarithm has a constant imaginary part. 

As explained in Ref. [|l6|], the zeros of the partition function should only appear outside the 
blue lines in Fig. ||. This can be inferred from the representation 

<j> dfi(dZ/dp)/Z = anYnJC) , (2.3) 
Jc Y 
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where n q (C) is the number of zeros of order q inside C. For large N, 

I dp(dZ/dp)/Z°c I dM 2 {dp/dM 2 )(M 2 -\/p) (2.4) 
Jc Jo 

The second term has a pole at j3 = 0, but it is compensated by a pole in M 2 because for small jS 
(large M 2 ), we have /3 ~ l/M 2 (see (2.2)). At infinite volume, the poles of (dfi/dM 2 ) are in the 
cut (the real interval [—8,0]). If the contour C' in the M 2 plane does not cross the cut, then there 
are no zeros of the partition function inside the corresponding C in the Z?-plane. We conclude that 
in the large-A^ limit, there are no Fisher's zero in the image of the cut M 2 plane. Fig. ^ illustrates a 
situation that we expect to be generic: the RG flows admit an analytical continuation in the complex 
plane until Fisher's zeros appear. 

We also considered the two lattice matching [|D|]. In short, we considered the sums of the 
spins in four L/2 x L/2 blocks B. NB denotes a nearest neighbor block of B. We define: 

((ZxeB$x)(LyeNB<i>y)) . 

R(P,L) = j _ ! P . (2.5) 

\{LxeB^x)(ZyeB^y))) 

A discrete RG transformation mapping j8 into j3' while the lattice spacing changes from a to 2a 
is obtained by matching: 7?(j8,L) = /?(/$', L/2). This method bypasses the difficult construction 
of the effective hamiltonian and the calculation of the field rescaling. The solutions were found 
numerically by Newton's method. When several solutions could be found, we picked the closest to 
the original /3 . We defined the ambiguity as | /3 — fS c losesi 1/10 — fhd.closest I • The results are presented 
in Fig. |forL = 4. 

It should noted that the appearance of a fixed point on the real axis near j8 ~ 0.3 is a finite size 
effect. As L increases, the fixed point moves to arbitrarily large /3 . In general, finite size effects in 
the complex plane are very intricate since the mapping of Eq. (^2|) becomes a ratio of polynomial 
and requires a multi-sheet Riemann surface to be inverted. This will be discussed at length in Ref. 

Ph. 
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Figure 3: Complex RG flows with the 2-lattice matching method for L = 4; Color scale :-Ln(ambiguity) 
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We also applied the two-lattice matching for Dyson' s hierarchical model with a Ising measure. 
It is a lattice model with block interactions depending on the details of the block configurations in 
a minimal way. In this model, the local potential approximation is exact. Its recursion formula 
is related to Wilson's approximate recursion formula (that allowed the first numerical RG calcu- 
lations) but the exponents are different. The probability distribution for the total spin in blocks 
can be calculated iteratively with this formula. The model has a continuous parameter that can be 
tuned in order to reproduce the scaling of a D-dimensional massless Gaussian field. For D = 2 
Dyson has proven rigorously the absence of transition [19]. For D = 3, it has a Wilson fixed point 
near j8 ~ 1.179. These facts are reviewed in Ref. j pO[ ] where the question of the improvement of 
the hierarchical approximation is also discussed. Complex RG flows showing bifurcations around 
lines of Fisher's zeros are shown in Fig. 2 of Ref. [|l3|]. The finite size effects and the ambiguity in 
the search will be discussed in Ref. [21 ]. 



3. Fisher's zeros in AD LGT 

At this point, we have not constructed complex RG flows for gauge theories, but we have 
designed numerically stable methods to calculate the Fisher's zeros in U(l) and SU (2) 4D lattice 
gauge theories. We use the spectral decomposition 

Srrmx „ 

dSn(S)e-P s (3.1) 



o 

with n(S) the density of states and jV the number of plaquettes. For large JV, we have 

n(S)e- pj " s = e-^VW-M = e ^{fW+(W(so){s-so) 2 +-) ( 3 _ 2 ) 

with s = S/-yV and /'(sq) = j8. f(s) can be interpreted as a color entropy density. If Ref" (so) < 0, 
the distribution becomes Gaussian in the infinite volume limit. Since Gaussian distributions have 
no complex zeros [|22|], the level curve Ref"(so) = is the boundary of the region where Fisher's 
zeros may appear. They typically delimit thin elongated regions ending at a complex zero of /". 



Results for SU (2) are shown in Fig. 3 of Ref. [13] and can be compared with the U(l) case in 



Fig. |J The main difference is that the imaginary part of the lowest zero appear to stabilize at a 
finite, non-zero value for SU (2) while this quantity goes to zero in the U(l) case. The details of 



the numerical construction of f(s) will be discussed in Ref. fl23Q. 

In the U (1) case, multicanonical methods were used [|4j] and naive histogram reweighting 
works well. The numerical error 8Z can be estimated from (n,-(5) — < n(S) >), where i is an 
index for independent runs. Zeros can be excluded if \SZ\ << \Z\. For SU (2), the imaginary part of 
Fisher's zeros are too large to use simple reweighting methods. By using Chebyshev interpolation 
for f(s) and monitoring the numerical stability of the integrals with the residue theorem as in Eq. 
(2.3), it is possible to obtain reasonably stable results. Unlike the U(\) case, the imaginary part of 
the lowest zeros does not decrease as the volume increases, but their linear density increases at a 
rate compatible with L~ 4 . The effect of an adjoint term (+0.5) is that the lowest zero goes down by 
about 40 percent. 



5 



Fisher's zeros and RG flows in complex coupling space 



Y. Meurice 



U(l) zeros 



L-4 f'-0 
L=4 res. 

L-6 f ' ' -0 
L=6 res . 
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Figure 4: Images of the zeros of f"(s) in the j3 plane (open symbols) and Fisher's zeros (filled symbols) for 
{/(!) on 4 4 (squares) and 6 4 (circles) lattices. 



4. Conclusions 

In summary, we have shown that it is possible to extend various RG flows to the complex j8 
plane. When the size of the system is comparable to the Compton wavelength of the gap, different 
methods can give very different answers. In all cases considered, the Fisher's zeros control the 
global behavior of the RG flows. Confinement occurs when the zeros leave an "open gate" on the 
real axis. We plan to study related questions for QED and SU (3) with various Nf. 
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